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Abstract 

Sum rules for the variation of finite-density spectral density of vector channel 
with baryon density are derived based on dispersion relations and the oper- 
ator product expansion. These sum rules may serve as constraints on the 
phenomenological models for the finite-density spectral densities used in the 
approaches motivated from QCD. Applying these sum rules to the rho meson 
in nuclear medium with a simple pole-plus-continuum ansatz for the spectral 
densities, we found that the qualitative features of the QCD sum-rule predic- 
tions for the spectral parameters are consistent with these sum rules; however, 

the quantitative QCD sum-rule results violate the sum rules to certain degree. 
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Study of hadronic properties at finite baryon density and temperature is essential in 
understanding the structure of matter. Recently, this subject has attracted much attention, 
motivated by the experimental attainment of hot and dense matter in heavy-ion collisions 
and by the theoretical expectations of a phase transition of matter from a hadronic phase 
to a quark-gluon plasma at high temperatures. 

Despite the difficulties due to the nonperturbative features of QCD at large distances, 
one may study the properties of hadrons and the QCD vacuum by investigating the two- 
point correlation functions of currents, carrying the quantum numbers of the system under 
study. This approach is based on the analytic properties of the two-point correlation func- 
tions and on asymptotic freedom. The hadronic spectral properties (e.g., masses, coupling 
constants, etc.) appearing in the spectral densities can be related via dispersion relations 
to the correlation functions evaluated in terms of quark and gluon degrees of freedom. In 
practical applications such as QCD sum-rule method [fl|,0 and analyses of lattice QCD data 
|| and interacting instanton model calculation 0], one needs to parametrize the spectral 
functions with a small number of parameters and to evaluate the correlation functions ap- 
proximately [e.g., operator product expansion (OPE), lattice simulations, and interacting 
instanton approximation] . The success of such approaches depends on correct understanding 
of the qualitative features of the spectral functions and accurate evaluation of the correlation 
functions from QCD. 

In this paper, we derive sum rules for the variation of finite-density spectral densities 
of vector channel with baryon density on the basis of dispersion relation and the OPE. 
These sum rules can be regarded as constraints on the phenomenological models for the 
finite-density spectral densities. We also apply these sum rules to the rho meson in nuclear 
medium with a simple pole-plus-continuum model for the spectral densities. We find that 
the qualitative features of the QCD sum-rule predictions for the spectral parameters are 
consistent with these constraints; the quantitative QCD sum-rule results, however, violate 
the constraints to certain degree. 
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Let us start with the correlation function of vector current at finite baryon density and 
zero temperature [0.0: 

U^(q,p B ) = % J d 4 xe^(TJ M (a;)J,(0)) Ps , (1) 

where J p (x) is a vector current constructed from light quark fields [e.g., J M (x) = 
| (u(x)'y l j i u(x) ± d(x)jpd(x)j]. Throughout this paper, the up and down quark masses are 
taken to be equal and all the quark fields in the current are assumed to have the same mass. 
The notation (• • -) PB denotes the expectation value on the finite-density ground state charac- 
terized by the baryon density p B in the rest frame and the four- velocity . In medium, there 
are in general two independent invariants in the vector channel, corresponding to the trans- 
verse and longitudinal polarizations Ht{q, Pb) and H(g, p B ). For simplicity, we will work in 
the rest frame of the medium, where = (1, 0), and take the three momentum to be zero 
q = 0. Then, since there is no specific spatial direction, the transverse polarization is related 
to the longitudinal one, Yl t (q ,p B ) = q^Hi(q , p B ), where Yli(q , p B ) = Il^(q , p B )/ (Sq^) 
f7]j5|||. The longitudinal part, IL(g ,p B ) = Hi(qQ,p B ), satisfies the standard dispersion 
relation @|| 

mips) - rds p -^4, ( 2 ) 

Jo s — qfi 

where p(s,p B ) = n^ImU^s, p B ) is the finite-density spectral density. Here we have omitted 
the subtraction terms, which can be eliminated by taking derivatives of both sides of Eq. (^) 
with respect to q 2 ,. The sum rules to be derived are, however, independent of this process. 

Using Eq. (Q), one can write the difference of the correlation functions evaluated at 
different baryon densities as 

An,(Q 2 ) ^ n,(Q 2 , PB ) - n,(Q 2 , P ' B ) = r ds P^pB)~ ^p'b) (3) 

Jo s + Q A 

where Q 2 = —q 2 and p' B denotes a different baryon density from p B . 

At very short distances, or at very high energies, the difference between two correlation 
functions with different densities should go to zero due to asymptotic freedom of QCD. We 
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can now look for the consequences of this statement for the difference of spectral densities. 
Expanding the right-hand side of Eq. @ for large values of Q 2 we get 

,2 



poo 

AU l (Q 2 ) = / ds[p(s,p B )-p(s,p' B )} 
Jo 



1 s s 



(4) 



On the other hand, for large Q 2 (i.e., in the deep Euclidean region), one can evaluate 
the correlation functions by expanding the product of currents according to the operator 
product expansion, which leads to 

n l (Q 2 lPB ) = Y,c n (Q 2 )(d n ) PB , (5) 

n 

where C n (Q 2 ) are the Wilson coefficients and O n are local composite operators constructed 
from quark and gluon fields. Here we have suppressed the dependence of both the coeffi- 
cients and the operators on the normalization point p. The operators O n are ordered by 
dimension (measured as a power of mass) and the C n (Q 2 ) for higher- dimensional operators 
fall off by corresponding powers of Q 2 . The Wilson coefficients only depend on QCD La- 
grangian parameters such as the quark masses and the strong coupling constant; all density 
dependence of the correlation function is included in the condensates (O n ) PB 0. Thus, one 
can express AUi(Q 2 ) as 

ATl l (Q 2 )=Y,C n (Q 2 )A(d n ) , (6) 

n 

where 

A(6 n ) = (O n ) PB - {d n ) p , B . (7) 

Note that the pure perturbative contribution [corresponding to the unit operator term, 
O n = 1, in the OPE] to the correlation function is independent of density, and thus does not 
appear in the difference ATli(Q 2 ). The lowest order contribution to AIl^Q 2 ) then comes 
from the condensates with lowest dimension (O n ^ 1), which, in the vector channel, are 
dimension four (including quark masses) condensates. Since AUi(Q 2 ) has dimension zero, 
the lowest order term in the OPE of AUi(Q 2 ) must be proportional to 1/Q A - However, the 



lowest order term in the phenomenological representation Eq. (£|) is proportional to 1/Q 2 - 
Therefore, we conclude that 



This is a rigorous result. In the OPE framework it simply follows from the observation that 
the pure perturbative contribution is density blind and the lowest-dimensional condensates 
have dimension four. Although the finite-density spectral densities cannot be measured 
directly from experiments, Eq. (§) will constrain the change of the spectral density with 
baryon density. Phenomenological parametrizations often used in applications such as QCD 
sum-rule calculations or analyses of lattice QCD data, must satisfy this constraint. If one 
adopts a pole-plus-continuum ansatz for the spectral density, Eq. (H) indicates that the 
change in the coupling for the pole is equal to the shift in the continuum. 
The OPE of Ui(Q 2 , pb) takes the general form 



where the ellipses denote the contributions of condensates with higher dimensions, the su- 
perscript indicates the dimension of operators, and the sum is over all contributing operators 
with a given dimension. The coefficients cf^ are dimensionless, and can, in principle, be 
dependent on Q 2 due to the perturbative corrections [only through strong coupling constant 



coupling constant and to the first order in quark masses, where the coefficients c\ become 
independent of Q 2 . 

We can then rewrite Eq. @ as 




(8) 





2 ) and m 2 /Q 2 ] [|I|. For simplicity, we will work only to the lowest order in the strong 




(10) 



Comparing the coefficients of 1/Q 4 in Eqs. (|J) and (0), one obtains 




(11) 
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Similarly, equating the coefficients of 1/Q e , one finds 

r s 2 [p(s, p B ) - p(s, p' B )\ ds = Y, of ) A(Of) . (f 2) 

Following the same pattern, one may derive an infinite series of sum rules, one for each 
OPE term (with fixed dimension) at small distances. The variation of the spectral density 
with baryon density must satisfy these sum rules. If one takes p' B = (i.e., in vacuum), 
these sum rules will constrain the change of the finite-density spectral density relative to 
the corresponding vacuum spectral density, which in some cases are experimentally accessi- 
ble. In the QCD sum-rule applications or the analyses of lattice simulation data, one often 
parametrizes the spectral densities with a pole representing the lowest resonance plus con- 
tinuum contribution roughly approximated by a perturbative evaluation of the correlation 
function, starting at an effective continuum threshold. One may apply the above sum rules 
to test this simple parametrization. 

In principle, the sum rules described here can also be used to determine the finite-density 
spectral density provided that its corresponding vacuum spectral density and the values of 
in-medium and vacuum condensates are known. In practice, however, one has to truncate 
the OPE as the number of condensates with the same dimension appearing in higher order 
terms become larger and there is no reliable way to estimate these condensates. As a result, 
one can only expect to utilize the first few sum rules (at best), which, along with a simple 
ansatz for the spectral density, may give rise to an estimate of the finite-density spectral 
parameters. 

We notice, however, that the sum rules of higher order (resulting from higher powers of 
1/Q 2 ) are more sensitive to the difference of the spectral density in higher energy region. 
In the QCD sum-rule approach, the Borel transformation suppresses the contribution from 
the higher energy region (continuum), though it introduces an auxiliary parameter (i.e., the 
Borel mass). Consequently, the spectral integral is saturated by the lowest resonance; the 
roughness of the approximation for the continuum is expected to have only small impact on 
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the spectral parameters for the lowest resonance The sum rules discussed here do not 
depend on any auxiliary parameter. However, for the sum rules to be useful in determining 
the spectral properties of hadrons, one needs to have a reliable model for the continuum. 

The sum rules derived in the present paper may look like the usual finite energy sum 
rules P,|9|JT0|1. However, we emphasize that in our study, we focuses on the difference of the 
correlation function evaluated at two different baryon densities, instead of the correlation 
function at a particular density. The present sum rules are for the variation of the finite- 
density spectral densities with baryon density, instead of the spectral density at a particular 
density. Our derivation relies on the subtraction procedure and on the asymptotic freedom 
of QCD, which allows a short distance expansion of both the correlation function and the 
phenomenological dispersion relation. This technique has been used by Kapusta and Shuryak 
|TT| in deriving the Weinberg-type sum rules at zero and finite temperature. One may derive 
analogous sum rules by using other type of subtraction scheme, instead of different baryon 



density (see for example Ref. [0). One can also extend the sum rules to other channels, 



which will be documented in Ref. [19 



We now turn to apply the sum rules Eqs. (g), and (|T|) to the rho meson in nuclear 
medium, where the vector current is given by J^(x) = | u(x)7 m m(x) — d^)^ ^d(x) fl]]. Since 
the complete spectral density in nuclear medium is not known experimentally, one cannot 
test directly whether the sum rules are indeed satisfied. 

Various investigators |||| have studied the properties of the rho meson in nuclear medium 
within QCD sum-rule approach. It is found that the rho meson mass (pole position of the 
propagator), the coupling of the vector current to the rho meson, and the effective continuum 
threshold all drop in nuclear medium. At nuclear matter saturation density, the rho meson 
mass decreases by ~ 15 — 18% relative to its vacuum value. One can test whether these 
in-medium spectral features predicted from QCD sum-rule calculations are consistent with 
the sum rules derived in the present paper. 

The explicit OPE result for the correlation function can be found in Refs. Pfl. To 
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the first order in the nucleon density pn, the in-medium condensates can be written as 
(0)pn = (^)vac + (O)nPn, where (O)n is the spin averaged nucleon matrix element |5|,[l4 



This linear approximation to the condensates is expected to be reasonable up to the nuclear 
matter saturation density |15|,|5|. Up to dimension six and to the linear order in p^, the 
result for the difference AIl^Q 2 ) can be written as 

77? I (y 

m{Q 2 ) = T1i(Q 2 ,Pn) - MQ 2 ,p'n = 0) = ^(qq)NPN + ^( -fG^GTUpN 

M 224 ^ A/f^ 

+^i A 2 +d PN ~ ^™ s (q q )Uw)NPN - ^ A T d pN , (13) 

where m q is the average of up and down quark masses, is the nucleon mass and A^ d {= 
A^ + A^) is a moment of the parton distributions in the deep inelastic scattering A^(p 2 ) = 
2 Jq 1 dxx n ~ l [q(x, p 2 ) + (— l)"g(x, /i 2 )], where q(x,p 2 ) and q(x, p 2 ) are the scale-dependent 
distribution functions for quarks and antiquarks (of flavor q) in the nucleon. Here the 
in-medium factorization approximation has been used for the four-quark condensates [i.e., 
(qq) 2 p N ~ (<?<?)vac — 2 vac (<?<?) nPn] 0- in Eq. (|13|), we only retained the terms included in 
Refs. |,| . 

Since the sum rules Eqs. (|j) and (|Tl~H r^) involve integrations of the spectral densities 



with different powers of s, it is likely important to incorporate the finite widths of the 
resonances and to have a reliable model for the continuum. Here we will adopt the pole- 
plus-continuum parametrization for both vacuum and in-medium spectral density, which is 
widely used in the QCD sum-rule calculations [ffl[5lJl6l|, and hence only expect to test the 
qualitative features of the finite-density spectral parameters. The detailed account of the 
finite width of the rho meson, along with the incorporation of a more reliable continuum 



model, will be given in Ref. |T3 



In the pole-plus-continuum approximation, one can write the vacuum spectral density 
as [0 

Pvac(s) = F5(s - m 2 ) + ^ (l + ^) 9(s - s ) , (14) 



and the finite-density spectral density as 



p(s, PN ) = Psc S(s) + F*5(s -rf) + A-(l + < ^L 



Is — s 



0) ' 



(15) 



where F (F*), m p (m*) and Sq (sq) are the coupling, the rho meson mass and the contin- 
uum threshold, respectively, and p sc denotes the contribution of Landau damping (or the 
scattering term) J7||| . 

Substituting the spectral ansatz Eqs. (p!4| -|T5|) and the OPE results into Eqs. (H) and 
(pTHT2|), we obtain the following sum rules 
1 



p sc + F* — F 



1 + — 







(16) 



m, 



Fmi 



16tt 2 



1 + — 

7T . 



(4 - sf) 



-m q (qq) N p N 
1 ,a. 



Ma 



A 2 p N 



m. 



Fml + 

p 24tt 2 



1 + — 

7T . 



(» 

224 
~8T 



- -LG^GnNPN- — AU+d 

Z4 7T 



na s (qq) vac (qq) N p N - ^^A^'Vat 



(17) 



(18) 



where we have assumed that the in-medium continuum threshold Sq is less than its vacuum 
value s . Substituting the QCD sum-rule predictions for the spectral parameters into these 
sum rules one may check how well these sum rules are satisfied. Alternatively, one may 
extract the in-medium spectral parameters by solving these three equations and compare 
the results with the QCD sum-rule predictions. Here we follow the latter. 

To obtain the finite-density spectral parameters m p *, F*, and Sq from Eqs. (|16|— 118 



one needs to know the various nucleon matrix elements appearing on the right hand sides 
as well as the corresponding vacuum spectral parameters. The nucleon matrix element 



(qq)N is related to the nucleon sigma term {qq)n = crjy/2m q ; we take <jn — 45MeV |L7| and 



5.5 MeV inyyl]. For the gluon matrix element, we use {(otg/^G^G^) N ~ —650 MeV 



|15|, 14] . The moments of parton distribution are taken to be A% +d ~ 0.938 and A u A +d ~ 0.121 



(at /i 2 = lGeV 2 ) 



We adopt (qq) y 



245 MeV) 3 [0 and a s ~ 0.3 § in our 



calculations. The nuclear matter saturation density is taken to be po = (110 Me V) 3 . We fix 
m p = 770 MeV and s = 1.5 GeV 2 [|IJ, and parametrize the scattering term as p sc = a p7v- 

In Fig. |I|, the resulting ratio of the in-medium rho meson mass to its vacuum value 
is plotted as function of nucleon density for different values of ao and a fixed F value 



F = 2f£ with f n = 93.5 MeV, which is obtained by using F = m 2 p /g 2 p gg] with the KSFR 
relation g 2 = m 2 /2f 2 . One notices that the in-medium rho meson mass drops relative to 
its vacuum value. Similar behavior is also found for both the coupling and the continuum 
threshold. These qualitative features are consistent with those predicted from the QCD 
sum-rule calculations. 

On the other hand, we observe that the quantitative result for the ratio is sensitive to 
the value of ao (i.e., p sc ). For ao = 0, the in- medium rho meson mass is only few percent 
smaller than its free space value even at the nuclear matter density. For ao = 3.6 GeV -1 , the 
rho meson mass drops to ~ 0.8m p at the saturation density. In Ref. the scattering term 
is neglected (i.e., ao = 0) while in Ref. || it is taken to be ao = l/2Mjv- The ratio m*/m p 
is found to be ~ 15 — 18% at p^ = po i n these references. However, we find that to obtain a 
15 — 18% decrease in the rho mass at the nuclear matter saturation density, one needs to use 
a value of ao ~ 3 GeV -1 , which is much larger than that used in Ref. ||. This inconsistency 
signals that the simple pole-plus-continuum model for the spectral densities with the QCD 
sum-rule predictions for the spectral parameters violates the constraints Eqs. (f|) and (|Tl|- 
12) to certain degree. A improved model beyond this simple parametrization for the spectral 



densities may be needed to satisfy the constraints. 

In Fig. ^j, we plot the ratio m*/m p as a function of the nucleon density for different 
values of F with fixed ao = 1/2Mn- It can be seen that the ratio is very sensitive to F, in 
particular for small F values. (For F > 0.035 GeV 2 , there is no real solution). Again, we 
note that a much smaller value of F is necessary to reproduce the QCD sum-rule result. It 
is also found that our result is relatively insensitive to the values of m p and so- 

In conclusion, we have derived sum rules for the variation of the finite-density spectral 
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density of vector channel with baryon density within the framework of operator product 
expansion and the dispersion relation. These sum rules may serve as constraints on the 
phenomenological models used in the QCD sum-rule calculations or in the interpretation 
of the lattice QCD data. We also noted that in principle one can use these sum rules to 
determine the qualitative properties of the finite-density spectral parameters if the corre- 
sponding vacuum spectral density and the in-medium and vacuum condensates are known. 
We applied the first three sum rules to the rho meson in nuclear medium with a simple 
pole-plus-continuum parametrization for the spectral densities, and found that the qualita- 
tive features of the QCD sum-rule predictions are consistent with our sum rules. However, 
the quantitative result shows that the simple ansatz with QCD sum-rule predictions for the 
spectral parameters violates the sum rules to certain degree. 

This suggests that the inclusion of the finite widths of the resonances and the refinement 
of the continuum model may be important. This point, along with the full detail of the 



present paper and its extension to other channels, will be reported elsewhere fL3||. 
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FIGURES 

FIG. 1. Ratio m*/m p , obtained from solving Eqs. (|i~6|-18), as a function of the nucleon density. 



The four curves correspond to ao = (solid), ao = 1.2 GeV 1 (long-dashed), ao = 2.4 GeV 1 
(dotted), and ao = 3.6 GeV -1 (dashed). The other input parameters are described in the text. 



FIG. 2. Ratio m*/m p , obtained from solving Eqs. (|16|4Tq), as a function of the nucleon den- 
sity. The four curves correspond to F = 0.005 GeV 2 (solid), F = 0.01 GeV 2 (long-dashed), 
F = 0.02 GeV 2 (dotted), and F = 0.03 GeV 2 (dashed). The other input parameters are the 
same as in Fig. H. 
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